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ABSTRACT 


A  computer  program  has  been  developed  which  follows  Light'  s 
phase- space  theory  of  three-body  ion- molecule  reactions.  The  program 
predicts  cross  sections  for  10  possible  reaction  branches  as  functions  of 
the  relative  kinetic  energy  of  the  colliding  particles  in  the  range  from 
thermal  energies  up  to  20  eV.  It  also  yields  partial  cross  sections  for 
the  populations  of  the  vibrational  levels  of  the  diatomic  molecule  in  the 
product  channel.  Included  in  the  results  are  reactions  for  the  following 
ion -molecule  pairs:  [O*  +  N2],  [N"*"  +  O2],  [O  +  N2  +  ]>  [He+  +  N2],  and 
[(0*k)+  +  O2].  The  effects  of  varying  the  level  of  excitation,  both  vibra- 
tionaily  and  electronically,  of  the  reactants  are  presented.  Increasing 
the  excitation  level  always  results  in  greater  cross  sections  for  reactions 
which  are  otherwise  endothermic  and  slightly  smaller  cross  sections  for 
reactions  which  are  otherwise  exothermic. 

All  endothermic  cross  sections  reach  maximum  values  as  functions 
of  the  barycentric  kinetic  energy  of  the  reactants.  All  exothermic  cross 
sections  decrease  with  increasing  barycentric  kinetic  energy  at  rates 
faster  than  predicted  by  the  Gioumousis  and  Stevenson  theory.  This  is 
due  to  the  greater  competition  between  the  products  for  the  available  phase 
space  which  occurs  at  elevated  energies.  Good  agreement  is  found  be¬ 
tween  the  theory  and  experiment  in  the  high  energy  region  above  5  eV  and 
in  some  cases  throughout  the  whole  energy  range.  Disagreement  is  most 
pronounced  for  low  energies.  It  is  felt  that  this  indicates  the  presence  of 
a  small  activation  energy  or,  in  the  case  where  the  products  differ  from 
the  reactants  by  a  charge  transfer  A"1"  +  BC  -*  A  +  BC+,  the  neglect  of  in¬ 
cluding  resonance  forces  in  the  calculation. 
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FIGURES 


1.  Production  of  NO  from  the  reaction  O  +  N2  •  The 
experimental  NO"1"  production  results  of  Giese^l  and 
of  Stebbings  _et_  al.^  are  shown  by  curves  XX  and 
YY  respectively, 

2.  Electronic  distribution  of  the  products  born  of  the 

reaction  0+  +  .  The  products  involved  in  dis¬ 

sociation  (D)  are  probably  2N  +  O"*".  The  theory  does 
not  always  predict  which  particle  is  ionized. 
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Electronic  distribution  of  the  products  from  the  re¬ 
action  (0+)5l<  (^D)  +  NO+.  The  increase  in  the  number 
of  products  follows  from  the  greater  energy  available 
when  one  of  the  reactants  is  in  an  excited  state.  More 
channels  are  exothermic  and  therefore  competitive  for 
tjhue  available  phase  space. 


Electronic  distribution  of  the  products  from  the  re¬ 
action  of  0+  +  N2  (v=4).  The  increased  vibrational 
excitation  of  the  reactant  molecule  produces  a  larger 
cross  section  (de- excitation)  for  formation  of  0+  +  N2 
(same  products  as  reactants)  because  the  de-excitation 
is  exothermic. 


5.  Vibrational  distribution  of  the  products  NO^"  (v=0,  1, 

.  .  .  ,  11)  +  N  from  channels  v=0,  1,2,3  are  exothermic. 
Channels  v  2  4  are  endothermic,  17 

6.  Vibrational  distribution  of  the  products  NO+  (v=0,  1, 

.  .  .  ,  11)  +  N  from  the  reaction  0+  +  N2  (v=4) .  Channels 
v  s  7  are  exothermic,  v  >  7  are  endothermic.  18 

7.  Products  NO+,  0+,  C>2+  from  the  reaction  N+  +  O^. 

The  experimental  results  shown  are  those  of  Stebbings 
et  al.^  Extrapolation  is  made  to  the  low-energy  experi¬ 
mental  data  of  Fehsenfeld  et  al.^®  19 
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Distribution  of  the  products  02*»  (O*®)*, 

and  from  the  reaction  +  O2.  The  ratios 

of  A1  to  A2,  A3  to  A4,  B1  to  B2 ,  and  B3  to  B4  are 
roughly  2  to  1 .  This  follows  from  the  phase- space 
theory  that  0^0*®  occupies  twice  the  phase  space  of 
O !®0 1® .  That  is,  reactions  of  the  type  A  +  B2  -*  AB 
+  B  have  cross  sections  2<T. 
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1.  INTRODUCTION 


Experimental  data  giving  cross  sections  as  functions  of  the 
reactant  ion  energy,  appearing  in  the  literature  on  ion-molecule  reactions, 
have  up  to  now  been  compared  only  with  the  theory  of  Gioumousis  and 
Stevenson  (G  and  S).^  This  theory  is  based  on  the  Langevin  model  which 
yields  orbits  of  a  charged  particle  moving  in  the  induced  dipole  potential 
which  exists  between  the  charged  and  neutral  particles.  When  the  impact 
parameter  is  less  than  a  certain  critical  value,  an  inward  spiraling  orbit 
takes  place.  Hence  if  the  impact  parameter  is  less  than  critical,  a  re¬ 
action  can  be  expected  to  occur,  whereas  impact  parameters  greater  than 
critical  cannot  give  rise  to  reaction.  This  reaction  cross  section  is  given 
by 

cr(Etr)  =  Tr(2e2a/Etr)^,  (1) 

where  e  is  the  charge  of  an  electron,  a  is  the  polarizability  of  the  neutral 
particle,  and  Etr  is  the  usual  kinetic  energy  of  the  reduced  mass  particle 
in  the  center  of  mass  coordinate  system. 

While  this  theory  has  been  successfully  compared  to  many  of  the 
experimental  data,  it  has  not  been  able  to  predict  such  features  as  the 
appearances  of  maxima  in  the  cross  section  versus  energy  curves  or  the 
rapid  fall  off  of  the  reaction  cross  sections  with  higher  energy,  nor  has 
it  been  able  to  deduce  what  fraction  of  the  products  are  in  excited  vibra¬ 
tional  and  electronic  states.  In  fact,  there  is  no  way  that  the  G  and  S 
theory  can  explain  which  of  the  products  are  most  likely  to  be  found  after 
reaction  or  the  fact  that  in  some  cases  no  reaction  occurs  at  all. 

This  is  not  surprising  when  one  realises  that  this  theory  gives  only 
an  upper  bound  to  the  reaction  cross  section.  Hence  one  should  not  expect 
to  find  any  "structure"  in  the  cross  section- energy  curve  predicted  by  it. 
Because  calculations  of  cross  sections  for  ion-molecule  reactions  are  prac¬ 
tically  intractable  quantum  mechanically  and  because  many  experimental 
data  indicate  something  other  than  that  predicted  by  Gioumousis  and  Steven¬ 
son,  we  feel  that  the  phase- space  theory  of  Light^  deserves  special  atten¬ 
tion®*  9.  To  study  this  theory  we  have  investigated  the  following  ion- 
molecule  reactions: 

0+  +  N_  -*  products, 

N+  +  o2 


products , 


products, 


These  reactions  are  of  some  importance  in  the  atmosphere.  Before  giving 
the  results  we  shall  first  review  the  Light  theory. 


2.  DESCRIPTION  OF  THE  LIGHT  THEORY 


Light1  s  phase- space  theory  is  based  upon  the  postulate  that: 

"The  probability  of  formation  of  any  given  product  in  a  1  strong  coupling1 
collision  is  proportional  to  the  ratio  of  the  phase  space  available  to  that 
product  divided  by  the  total  phase  space  available  with  conservation  of 
energy  and  total  angular  momentum."  The  strong  coupling  is  necessary 
to  ensure  that  the  complex  formed  by  the  reactants  "loses  all  memory" 
of  its  initial  state,  and  consequently  that  its  decomposition  be  governed 
by  phase  space  available. 

The  classical  three-particle  phase- space  element  is  the  usual 

3  3  3 

dr  =  .n  d  r.d  p.  ,  (3) 

i=  1  l  l 

where  r.  and  p.  are  the  position  and  momentum  coordinates  of  particle  i. 
By  separating  out  the  center  of  mass  motion,  writing  JT  =  JoRB  +  JROT> 
and  taking  J -p  to  lie  along  the  z-axis,  we  find  after  a  Jacobian  transforma¬ 
tion,  the  reduced  phase- space  element 

dF  ■  <JT2  +  JROT2  -  2JTJZROT>'4  JTdJTdJZTdV“0RB 

dBTdtTdJROTdJZROTdaROTdf3ROTdEVIBdtVIB  ,  (4) 

where  the  symbol  J  refers  to  angular  momentum,  E  refers  to  energy,  t 
to  time,  a  to  the  angle  conjugate  to  J,  and  (3  to  the  angle  conjugate  to  the 
Z  component  of  J.  The  subscript  T  means  total,  ORB  means  orbital  and 
refers  to  the  orbital  motion  of  the  single  particle  about  the  molecule,  ROT 
means  rotational  and  refers  to  the  rotation  of  the  diatomic  molecule  about 
its  mass  center,  and  VIB  means  vibrational  and  refers  to  vibration  of  the 
diatomic  molecule.  ^  By  noticing  that  the  integrals  over  C*ORB’  ^ROT* 
PrOT»  Pt*  tT>  and  ^VIB  Yield  constants  independent  of  channel  and  vibra¬ 
tional  energy,  we  can  write  for  the  relevant  phase  space  available  with 
conservation  of_total  energy  E'p,  vibrational  energy  Eyjg,  and  total  angu¬ 
lar  momentum  J,j,, 

r  =  JJ  }  +  (JROt/Jt)  "  2JZROT//Jt]  ^  dJROTdJZROT  *  ^ 
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It  should  also  be  noted  that  the  integral  given  by  Eq.  (9)  of  Reference  7 
is  incorrect,  as  suggested  by  Pechukas  and  Light  in  a  later  paper. ^  The 
correct  expression  is  given  by  our  Eq.  (5). 

The  limits  of  integration  of  Eq.  (5)  are  determined  by  the  following 
physical  constraints: 

a.  Since  energy  must  be  conserved  we  have 


ET  =  E,r 


ROT 


ROT 


2 

where  E^r  is  the  translational  kinetic  energy  p,v  /2,  p  is  the 
reduced  mass  of  the  diatomic  molecule  in  channel  k,  and 

k 

Erot  *-s  rotational  energy.  This  is  given  by 


where  I  is  the  moment  of  inertia  of  the  diatomic  molecule. 
The  superscripts  o  and  i  refer  to  the  input  and  ith  output 
channel  respectively,  and  Q0j  is  the  exothermicity  for  the 
reaction.  Thus  an  endothermic  channel  will  have  a  negative 
Q.  Because  Etr  must  be  positive,  Eq.  (6)  yields  for  each 
channel  i  an  upper  bound  on  the  rotational  angular  momentum 

JROT2/21  S  ET  '  V  +  Qoi  =  C  •  (8) 


b.  Since  the  molecule  in  the  output  channel  must  be  stable,  we 
must  also  limit  the  molecule  from  dissociating  rotationally . 
Therefore 


jrot^1<d/  ■  <9> 

where  DV1  is  the  dissociation  energy  of  the  diatomic  molecule 
from  the  vibrational  level  v. 


c.  The  final  requirement  is  that  the  products  must  separate. 

When  the  products  are  an  ion  and  a  neutral  particle,  the  effec¬ 
tive  one- dimensional  central-field  potential  is  given  by 


V{r)  =  -  &  “ cl/ 2r^  +  E^b^/r2 
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where  a  is  the  polarizability  o£  the  neutral  particle,  Ej  is 
the  final  translational  kinetic  energy,  and  bf  is  the  impact 
parameter  for  the  products.  This  interaction  potential  can  be 
expected  to  be  physically  valid  only  for  long  ranges.  Clearly 
the  short-range  repulsion  is  not  included;  however,  it  would 
only  be  important  for  large  values  of  the  mass  reduced  orbital 
angular  mom«»tum,  JORB? =  Efbf  .  For  low  relative 
kinetic  energies,  the  centrifugal  barrier  occurs  in  the  region 
of  the  attractive  induced  dipole  potential.  For  intermediate 
values  of  E£,  quantum  mechanical  resonance  forces  must  be 
taken  into  account.^  Present  calculations  are  underway  which 
include  resonance  potentials.  The  results  will  be  reported 
soon  for  several  pertinent  ion-molecule  reactions. 

The  charge- quadrupole  interaction  which  leads  to  a  potential 
varying  as  r"3  can  be  neglected  on  the  basis  that  it  is  usually 
much  smaller  than  the  induced  dipole- charge  interaction.  This 
has  been  investigated  by  Arthurs  and  Dalgarno.^  Applying 
their  calculation  to  H£  +  fT?”4  leads  to  the  cross  section 


°AD  =  1'08  (Etr) 


-2/3  <i2) 


w1  re  Etr  is  the  relative  kinetic  energy  of  the  reactants  in 


The  G  and  S  cross  section  for  this  reaction  is 


"GS  =  15-08  Kr)'4  (A2»  i 


hence,  it  is  clear  that  the  charge- quadrupole  interaction  is 
only  important  at  extremely  low  energies,  well  below  room 
temperature . 

The  restriction  that  the  products  must  separate  can  be  written 
very  simply  as 

Ef>V(r*)  ,  (13) 


where 


* =  (e2a/Efbf2)4  - 
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and 


❖ 

V(r  )  = 


Efbf  '  /(2e2a) 


(15) 


are  the  position  and  value  of  the  potential,  respectively,  for 
the  centrifugal  barrier.  This  can  be  converted  into  a  condition 
which  limits  the  range  of  values  taken  by 

Using  conditions  a  and  c,  we  find  for  the  total  phase  space  available 
to  products  i  in  vibrational  state  v,  the  result 

r(Jr  J2>  JT>  0  =  Ix  +  I2  S<-  +  I3  sgn(JT  -  j2) 

+  I4  sgn(JT  -  J{)  +  I5  S(J2  -  JT)  S(JT  -  (16) 

where  in  the  above 


S(A  -  B)  =  < 


1,  A  >  B 
0,  A  <  B 


sgn(A  -  B)  = 


+  1,  A  >  B 
-1,  A  <  B  , 


(17) 


and 


1 1  =  (8(1  2e2a) 


1/4 


2 

/  (- 


JROT  /2I) 


1/4 


dJ 


ROT  ’ 


h  ~  JT*JJ  +  2  J1  * 


I3  JTJ2  +  2  J2 

I4-JTIJ1|-5J12 


I-  =  -  J, 


(18) 


6 


The  values  Jj  and  3^  correspond  to  the  roots  Xj  and  X2  given 
in  Reference  7.  Briefly,  they  are  bounds  on  the  values  that  JroT  can 
assume  when  the  inequality  E£  a  V(r*)  is  applied.  They  are  given  by 

J  =  (210*  X.  ,  i  =  1,  2  ,  (19) 

1  1 

and  always  satisfy  the  inequality  J2  s  |j]J.  The  total  phase  space  available 
to  all  products  is  given  by  summing  over  both  vibrational  and  product 
(electronic)  states. 

r(ET,  jt)  sssr.dj,  J2,  JT,  0  .  (20) 

By  not  restricting  T  through  use  of  condition  b,  we  are  including  the  phase 
space  available  to  rotational  dissociation. 

Since  the  product  channel  contains  a  molecule,  the  product  >hase 
space  available  must  be  restricted  by  condition  b.  Writing  J3  =  (2IDV^^ 
for  brevity,  we  find  for  the  product  phase  space 

r.'  (jv  J2,  JT,  e)  =  r.Uj,  j2,  jt,  os(j3-j2) 

+  ri(Jr  J3*  jt’  0  S(J2-J3)  s(J3-|J1l) 

+  [(2jtj3-jt2)  s(j3-jt)  +  j32  s(jt-j3)] 

X  S  (-Jx)  SdjJ-^)  .  (21) 

The  probability  for  the  production  of  molecule  i  in  the  vth  vibra¬ 
tional  state  is  then  given  by 

p(et-  jt-  ev)  =  ri' /r(Er  jt)  •  <22> 

The  cross  section  for  this  production  is  then 

b 

max 

V  (Etr°,  i»  vj  =  2tt  J  P^Et,  JT,  EvX)  bdb  ,  (23) 

o 
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where  b  is  the  impact  parameter  in  the  reactant  channel  and 


b 

max 


2e2a/E 


1/4 


tr 


By  using  Eq.  (23),  we  are  assuming  that  the  initial  energy  and  angular 
momentum  of  rotation  of  the  reactant  molecule  is  negligible  in  comparison 
with  the  total  energy  and  total  jmgular  momentum  respectively.  This  al¬ 
lows  us  to  take  =  2^iEtr°b  .  We  point  out  that  if  P  =  1,  the  cross 

section  reduces  to  that  predicted  by  the  G  and  S  theory;  therefore,  the 
phase- space  theory  enables  us  to  predict  what  fraction  of  the  total  reaction 
cross  section  results  in  the  production  of  a  particular  product  channel  i 
and  vibrational  state  v. 


The  cross  section  for  dissociation  is  found  by  subtracting  the  sum 
of  the  formation  cross  sections  from  the  maximum,  or 


^DISS 


=  trb 


2 


max 


Z  Z  a(E,  °,  i,  v)  . 

i  v  '  tr  ' 


(24) 


Our  calculations  have  been  carried  out  on  a  fast  computer,  using  an 
expanded  version  of  Light'  s  original  program  to  include  10  output  channels. 
This  includes  a  sub-program  which  yields  the  roots  Xj  and  X2  of  Eq.  (19) 
in  Reference  7  and  calculates  the  integral  Ij  by  a  simple  numerical  quad¬ 
rature  . 


Our  calculations  make  use  of  many  data  already  available  in  the 
literature  on  atomic  and  molecular  constants.  Molecular  potential  energy 
curves,  interatomic  separations,  dissociation  energies,  and  vibrational 
energy  levels,  when  given, were  taken  from  the  work  of  Gilmore.  Atomic 
polarizabilities  were  taken  from  the  work  of  Parkinson,*®  and  molecular 
polarizabilities  were  obtained  from  the  experimental  results  of  Rothe  and 
Bernstein*^.  In  all  cases,  vibrational  levels  were  obtaihea  using  the 
formula*® 


Evk  =  (v  +  1/2)  h|/k  -  (v  +  1/2)2  (hl^)2/*  Dek  ,  (25) 

where  De  is  the  dissociation  energy  from  the  bottom  of  the  potential 
well  of  the  product  molecule  k  and  fitting  the  frequencies  to  the  poten¬ 
tial  energy  curves  given  by  Gilmore. 


8 


3.  DISCUSSION  OF  THE  RESULTS 


The  results  of  the  calculations  are  shown  in  Figures  1  through  19. 

We  have  plotted  cross  sections  in  units  of  10“  ^  Cm^  =  1  versus  the 
incoming  relative  kinetic  energy  of  the  reactants,  (jlv^/2,  in  eV  units.  For 
most  of  the  energy  range  shown,  a  simple  conversion  from  barycentric  to 
laboratory  ion  energies  is  found  by  writing 

Elab  =  (m/H°  E  ’  (26) 

where  m  is  the  incoming  ion'  s  mass  and  p  is  the  reduced  mass  of  the 
reactants.  For  example,  in  Figure  1  the  incoming  reduced  mass  is  p  = 
16*28/44  amu,  while  the  ion'  s  mass  is  16.  Hence  Eja^  =  1.57  E. 

The  straight  negative  sloped  line  labeled  G  and  S  in  all  of  the 
figures  represents  the  Gioumousis  and  Stevenson  cross  section  for  each 
set  of  reactants.  Thus  it  appears  as  shown  when  plotted  on  the  log- log 
graph.  Furthermore,  this  cross  section  is  the  maximum  obtainable  when 
the  induced  dipole  potential  is  used.  By  summing  over  all  production 
cross  sections,  one  obtains  the  G  and  S  result. 

Each  curve,  unless  otherwise  indicated,  represents  a  composite  in 
that  it  is  a  sum  over  all  possible  vibrational  levels  of  the  diatomic  molecule 
in  the  output  channel  depicted.  ^  Thus  Curve  A1  in  Figure  2  represents  the 
total  cross  section  for  NO*(X  E*)  production.  Figure  5  shows  the  vibra¬ 
tional  distribution  of  the  product  A1 . 

Curves  A2  and  A3  of  Figure  2  represent,  respectively,  the  produc¬ 
tion  of  the  excited  neutral  states  N(^D)  and  N(^P)  together  with  NO*(X  *E+) . 
Since  these  reactions  are  endothermic,  their  cross  sections  are  zero  at  low 
energies.  Curve  A  of  Figure  1  is  the  sum  of  Al,  A2,  and  A3,  and  thus  re¬ 
presents  the  total  N  +  NO*  production. 

Curve  B  of  Figure  1  represents  the  cross  section  for  production  of 
the  same  products  as  the  reactants.  For  low  energies,  this  cross  section 
is  non- reactive  in  that  the  products  are  also  in  the  same  vibrational  state 
as  the  reactants.  For  higher  energies,  however,  higher  vibrational  states 
will  appear.  This  effect  can  be  seen  in  Figure  5,  which  shows  the  vibra¬ 
tional  distribution  of  the  cross  section  for  the  production  of  NO*(X  *E* 
v=0,  1,  ...,  11).  In  all  the  figures,  at  higher  energies  all  vibrational 
levels  are  excited  and  the  cross  sections  are  seen  to  approach  each  other 
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asymptotically.  Low-energy  cross  sections  for  particular  vibrational 
levels  only  appear  for  those  vibrational  levels  which  are  exothermic. 

Comparing  Figures  7  and  10,  we  see  that  for  low  energies,  agree¬ 
ment  is  best  for  the  reaction 

N+(3P)  +  Oz(X  3S  ‘  v=0)  -  0(3P)  +  NO+(X  XE+)  , 

8 

while  for^energies  above  5  eV,  the  experimental  results  agree  with  the 
total  NO  production  (the  sum  of  Curves  A1  through  A6,  represented  as 
A  on  Fig.  7).  We  point  out  that  the  above  reaction  is  more  than  6-eV 
exothermic,  and  hence  one  may  not  expect  to  find  any  important  effects, 
due  to  activation  energies.  The  other  reactions,  A2  through  A6  in  Figure 
10,  are  not  as  exothermic  (A5  and  A6  are  endothermic);  therefore,  if 
activation  is  present,  these  reactions  may  not  be  observed  at  low  energies. 
This  could  explain  the  good  agreement  of  experimental  results  with  A1  at 
low  energies.  At  higher  energies,  the  experimental  results  best  fit  the 
total  NO+  production  shown  as  A  in  Figure  7.  This  agreement  is  reason¬ 
able  in  light  of  our  previous  discussion  of  phase- space  competition  at 
higher  energies. 

Discrepancies  loom  large,  as  a  usual  rule,  whenever  a  resonant 
charge-transfer  process  is  observed  to  take  place;  i.e.,  our  results  yield 
cross  sections  for  rearrangement  resulting  in  charge  transfer  which  is 
usually  smaller  than  observed  experimentally.  By  taking  into  account 
resonance  forces,  we  hope  to  be  able  to  predict  cross  sections  for  low- 
energy  charge  transfer  which  give  better  agreement  with  experiment. 

This  will  appear  in  a  forthcoming  paper.  In  the  event  that  the  reactant 
channel  "sees"  a  resonance  force,  we  could  expect  to  find  total  cross 
sections  which  would  exceed  the  G  and  S  maximum.^  In  the  event  of 
actual  rearrangement,  the  results  seem  promising.  Agreement  is  expected 
to  be  best  when  the  output  channels  are  highly  exothermic,  i.e. ,  when  the 
effects  of  any  activation  energy  are  expected  to  be  small.  Thus,  see  the 
results  for  the  reaction  N+  +  O2  -*  O  +  NO+.  Here  (Fig.  7)  good  agreement 
is  found  between  theory  and  experiment  throughout  the  energy  range  shown. 

The  most  significant  feature  of  the  theory  is  the  fact  that  our  re¬ 
sults  yield  cross  sections  for  exothermic  reactions  which  fall  off  faster 
with  increasing  barycentric  kinetic  energy  than  that  predicted  by  Gioumousis 
and  Stevenson  when  the  energy  exceeds  5  eV  or  so.  This  is  in  good  agree¬ 
ment  with  experimental  measurements^’  ^  ^  *  (witness  Figs  .  1,  7,  10, 

and  16).  Our  results  indicate  that  this  is  due  to  phase- space  competition. 

As  the  incoming  energy  increases,  more  channels  open  for  the  reaction 
path  to  follow.  Thus  each  channel  shares  in  the  output.  Finally,  when  the 
incoming  energy  is  high  enough,  the  exothermicities  become  small  portions 
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of  the  total  available  energy.  Hence  all  products  become,  roughly,  equally 
likely.  The  small  value  for  each  cross  section  is  attributed  to  the  fact 
that  more  phase  space  is  available  to  dissociated  products  at  high  energies 
than  to  stable  vibrating  molecules,  since  statistically  a  continuous  energy 
distribution  wins  out  over  a  discrete  one. 

The  disagreement  at  low  energies  is  significant  and  indicates  either 
that  an  activation  energy  and/or  "selection  rule"  restriction  is  in  operation. 
Bates  and  Lynn  point  out^^  that  even  though  a  reaction  may  be  favorable 
energetically,  it  may  be  ruled  out  if  an  electronic  transition  is  required  to 
"make  it  go".  Thus  they  would  predict  low  cross  sections  for  asymmetric 
resonant  systems.  In  the  reaction  O*  +  N2  “*  NO*  +  N  (Fig.  1),  this  is, 
however,  not  the  case.  The  dissociated  products  and  reactants  match, 
i.  e.  ,  O*  +  N2  -*  O*  +  N  +  N  *-  NO*  +  N,  hence  no  electronic  transition  is 
required.  Thus  we  are  at  a  loss  to  explain  the  experimental  data  on  this 
basis.  Recently,  Schmeltekopf  and  his  associates  at  the  National  Bureau 
of  Standards,  Boulder,  Colorado,  ^  observed  an  increase  of  about  a  factor 
of  20  in  the  reaction  rate  for  O*  +  N2  ,-^NO*  +  N  when  the  N2  was  initially 
in  a  higher  vibrational  state.  Our  results  are  shown  in  Figures  4  and  6. 
More  vibrational  levels  are  excited  in  the  NO*  channel  as  a  result  of  in¬ 
creasing  the  reactant  vibrational  energy.  The  measured  results  in  this 
case  do  seem  to  agree  with  that  predicted  here,  thus  confirming  the  specu¬ 
lation  that  the  theory  best  fits  reactions  which  are  strongly  exothermic. 

The  results  predicted  by  the  reaction  (Figs.  16-18) 

He+  +  N2  -  He  +  (n 2+)*  -  He  +  N  +  N+  (2  7) 

deserve  special  mention.  Here  there  is  enough  ionization  energy  available 
so  that  the  dissociative  reaction  appears  as  an  exothermic  channel  with  the 
cross  section  increasing  with  decreasing  energy.  This  is  in  contrast  to 
any  of  the  other  cases  shown,  where  dissociation,  being  an  endothermic 
channel,  does  not  take  place  until  the  translational  kinetic  energy  exceeds 
the  energy  defect.  Our  results  for  dissociative  charge  transfer  in  this 
case  fall  much  below  the  experimental  results  because  the  process  is  one 
in  which  resonance  forces  are  probably  quite  large.  Thus  the  accidental 
resonance  reaction 

He+(2S)  +  N2(X  1Sg+  v=0)  -  He^S)  +  N£+(C  v=4)  +  O  eV  (28) 

could  predominate,  as  may  be  deduced  from  the  theory  of  Rapp  and 
Francis.  ^  The  phase- space  theory  as  presented  will  always  predict 
cross  sections  smaller  than  the  G  and  S  theory,  and  such  a  reaction  as 
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(28)  constitutes  statistically  a  small  part  of  the  total  Gioumousis  and 
Stevenson  cross  section.  It  is,  however,  of  interest  to  see  the  vibrational 
distribution  of  the  output  cross  section  for  this  reaction,  and  this  is  shown 


in  Figure  18.  By  including  resonance  forces,  we  might  expect  the  cross 
section  for  production  of  the  v=4  level  of  N2+  C  to  be  much  larger 

than  shown  here.  The  potential  energy  curve  crossing^  of  the  ^11  state 
with  the  C  state  would  then  lead  to  dissociation  of  the  ^  +  molecule^, 


We  have  also  calculated  cross  sections  for  the  reactions  N0+  +  O 
and  (Oao)  +  O2 .  The  former  are  shown  in  Figures  13  through  15  and  the 
latter  in  Figure  19.  The  latter  reaction  is  of  interest  for  future  experi¬ 
mentation,  since  it  predicts  a  rather  simple  isotope  effect.  If  we  consider 
reactions  of  the  type  A+  +  B 2,  we  have  as  possible  products 


(a)  A  +  B2  + 

(b)  AB  +  B+  (2  ways  to  form  this) 

(c)  AB+  +  B  (2  ways  to  form  this) 

(d)  A+  +  B2 


Consider  products  (a)  and  (b) .  If  A  is  an  isotope  of  B,  the  reaction 
for  (b)  will  occur  statistically  twice  as  many  times  as  that  for  (a).  Hence 
rearrangement  charge  transfer  (b)  for  the  reaction  (0^)+  +  O2  -*0  4-  (O^C)18)  + 
will  occupy  twice  the  phase  space  available  to  simple  charge  transfer  (a). 
Thus  this  reaction  should  prove  a  good  test  for  the  phase- space  theory. 
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cr  (A2) 


ENERGY  (eV) 

REACTANTS  0  +  (4S°)  +  Nz  (X  1 25  v  =0) 


-f 

Fig.  l--Froduction  of  NO  from  the  reaction  O  +  N2.  The  experimen¬ 
tal  NO"*"  production  results  of  Giese^  and  of  Stebbings  et  al.  ^ 
are  shown  by  curves  XX  and  YY  respectively* 
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ENERGY  (eV) 

REACTANTS  0+(4S°)+N2  (X'  v  =0) 


4. 

Fig.  2- -Electronic  distribution  of  the  products  born  of  the  reaction  O  + 
N2.  The  products  involved  in  dissociation  (D)  are  probably  2N  + 
0+.  The  theory  does  not  always  predict  which  particle  is  ionized. 
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A2  N0+(X'2$)+N  (ZD) 
A3  N0*(X'2g)  +  N(3P) 
6  I  0+(4S)  +  N,(X'2:S) 
B2  0+(2D)  H-N2(X'2+ ) 

Cl  0(3P)+N|(Xz2$) 
C2  0(3P)  +  NJ(A27Tu) 
C3  oCdj+nkx2^) 
E  N*(3P)  +  NO(Xz77J.) 
D  DISSOCIATION 


ENERGY  (eV) 

REACTANTS  0+(2D)*+  N2  (X 1  2$  v  * 0) 


+ 

Electronic  distribution  of  the  products  from  the  reaction  (O  ) 
(^D)  +  NO*.  The  increase  in  the  number  of  products  follows 
from  the  greater  energy  available  when  one  of  the  reactants  is 
in  an  excited  state.  More  channels  are  exothermic  and  therefoi 
competitive  for  the  available  phase  space. 


ENERGY  (eV) 

REACTANTS  0+(4S°)  +  N2(X'  2$  v  =4) 

Fig.  4- -Electronic  distribution  of  the  products  from  the  reaction  of 

0+  +  N2  (v=4).  The  increased  vibrational  excitation  of  the  re¬ 
actant  molecule  produces  a  larger  cross  section  (de- excitation) 
for  formation  of  0+  +  (same  products  as  reactants)  because 
the  de- excitation  is  exothermic. 
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ENERGY  (eV) 


0+  (4S)  +NZ  (X1 2g  v  =  0)  — * 

N  (4S°)  +  N0+(X‘  X*v  =  0, 1,-, 1 1) 


Fig.  5-- Vibrational  distribution  of  the  products  NO+  (v=0,  1,  .  .  .  ,  11) 

+  N  from  channels  v=0,  1,  2,  3  are  exothermic.  Channels  v  >  4 
are  endothermic. 


ENERGY  (eV) 

0+(4S)  +N2  (X'Sg  v  =4)  — * 
N  (4S)+N0+(Xi  2*  v®0,l;",M) 


Fig.  6 --Vibrational  distribution  of  the  products  NO  (v=0,  1,  .  ..,  11) 
+  N  from  the  reaction  0+  +  N^  (v=4).  Channels  v  £  7  are  exo¬ 
thermic,  v  >  7  are  endothermic. 
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PRODUCTS  N.  \ 

A  0+N0+ 

EXP  0+N0+  \ 

D  DISSOCIATION 
E  0+(4S)+N0(X27Tr) 
B  N+(3P)4-02(X32g) 
C  N+0|(X277^) 

C  I  N  (4S)+0| 

C2  N  (2D)  +  0£ 


EXP 


ENERGY  (eV) 

REACTANTS  N*(3P)  +  02  v  =  0) 


4"  *4’ 

-Products  NO  ,  O  ,  from  the  reaction  N  +  O2 .  The  experi¬ 

mental  results  shown  are  those  of  Stebbings  et  al.  2  Extrapola¬ 
tion  is  made  to  the  low-energy  experimental  data  of  Fehsenfeld 
et  al. 


.01  0.1  I  10  100 


ENERGY  (eV) 

REACTANTS  N+( ' D)% 02 (X3 Zg  v  =0) 

Fig.  8--Product  distribution  NO+,  N+,  and  0+  from  the  reaction 

(N*)  ‘  +  C>2.  The  de- excitation  cross  section  is  given  by  B. 


20 


ENERGY  (eV) 

REACTANTS  N*(3P)  +  02  (X3Z5  v  *  4) 


Fig.  9 ■  -  Product  distribution  NO+,  N+,  02  +  ,  and  0+  from  the  reaction 
N+  4-  Oz  (v=4) 
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PRODUCTS 

A  I  0(*P)+  N0+(X'2+) 
A 2  0('D)+ NO  +  (X'E+) 
A3  0(3P)+ N0*(d32+) 
A4  0('S)+  NO  +  (X'2+) 
A5  0('D)+N0  +  (a32+) 
A6  0(3P)+0N  +  (  3A  ) 
EXP  0+N0+ 


ENERGY  (eV) 

REACTANTS  N+(3P)  +  0z  (X3Zg 


Fig.  10-- Electronic  distribution  of  the  products  O  +  NO  from  the  re¬ 
action  N+  +  O2 .  The  experimental  results  are  those  of  Stebbing 
aL  2  at  higher  energies.  These  have  been  extrapolated  to  the 
room  temperature  data  of  Fehsenfeld  et  al.^® 
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PRODUCTS 

A  I  0(3P)  +  NO+(X'2+) 
A2  0('D)+ NOTX'S*) 
A3  0(3P)+ N0+(d32*) 
A4  0('S)+ N0  +  (X'2+) 
A5  0('D)+N0  +  (a32+) 
A6  0(3P)+  ON  +  (  3A  ) 


ENERGY  (eV) 

REACTANTS  N+( 'D)*+ 02(X3  2g  v -0) 


Fig.  1 1-- Electronic  distribution  of  the  products  O  +  NO+  from  the 
action  (hr1")5'  +  O-. 


Fig.  12--Electronic  distribution  of  the  products  O  +  NO+  from  the  re¬ 
action.  N+  +  O2  (v=4) 
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ENERGY  (eV) 

REACTANTS  0  (3P)  + Nj  (X8S5  v«0) 


Product  distribution  NO  ,  O 


N- 


+ 


and  from  the  reaction 


1 — !~"l  I  I  I 


T  I  "I  I 


100 


PRODUCTS  A  I  N0+(X'2g)+N  (4S) 
A2  N0+(X'2J)+N  (*D) 
A3  N0*(X'5:J)+N  (3P) 
:  B  I  OVS)  +  N2(X'S5) 

B2  0+(2D)  +  NZ(X'S|) 

Cl  0  (3P)  +  N|(X22+) 
C2  0  (sP)  +  N2(A27Tu) 
C3  0('D)  +  N5(X225) 
Aj\  E  N+(3P)  +  NO(X277:) 
\p  DISSOCIATION 
Bl 


A2 


A3 


C2 


XI 


/ 


J _ I  I  I  I 


B2 


j _ i  — . j  i. 


"C3 


JJLL 


1,111 


0.1  I  10 

ENERGY  (eV) 

REACTANTS  0  ( 1  D)*+  N^X2^  v*0) 


too 


-}-  -j"  *f*  •{- 

Fig.  14--Product  distribution  NO  ,  O  ,  N2  ,  and  N  from  the  reaction 

o'  +  V 
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cr  (A2) 


ENERGY  (eV) 

REACTANTS  0  (3  P)  +  N|(XZ  2$  v  =4) 


•f*  "f"  *f~  “j~ 

Fig.  15 --Product  distribution  NO  ,  O  ,  ,  and  N  from  the  reaction 

O  +  N2+  (v=4) 
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<r  (A2) 


ENERGY  (eV) 

REACTANTS  Hj(2 S)+  N2  (Xz2$  v  =  0) 


Fig. 


+  + 

l6--Product  distribution  N£  ,  N  ,  and  He 
4-  N2-  The  experimental  data  for  N2  + 
from  Stebbings  et  al.^ 


from  the  reaction  He 
production  are  taken 
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.01  0.1  I  10 

ENERGY  (eV) 

He(2S)  +Nz(Xz2gv=0)— > 
He  ('S)+N|(C22Uv=0,I,*"7) 


Fig.  18- -Vibrational  distribution  of  the  products  (v=0,  1, 

+  He  from  the  reaction  He*  +  N2 


rrrrra 


I-J.I1.HJ 

100 


....  7) 


30 


ENERGY  (eV) 

REACTANTS  (0,8)+(4S°)+02(X3  Eg  v*0) 


Fig.  19- -Distribution  of  the  products  0£  + ,  (O  O  )  ,  (O  )  ,  and 

(018)+  from  the  reaction  (O^8)"1"  +  C^.  The  ratios  of  A1  to  A2 , 
A3  to  A4,  B1  to  B2,  and  B3  to  B4  are  roughly  2  to  1 .  This  fol¬ 
lows  from  the  phase- space  theory  that  O^O^8  occupies  twice 
the  phase  space  of  oV°Q^8.  That  is,  reactions  of  the  type  A  + 
B2  AB  +  B  have  cross  sections  2a, 
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